Magneto-thermodynamics of the spin-^ Kagome antiferromagnet 
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In this paper, we use a new hybrid method to compute 
the thermodynamic behavior of the spin-| Kagome antifer- 
romagnet under the influence of a large external magnetic 
field. We find a T 2 low-temperature behavior and a very low 
sensitivity of the specific heat to a strong external magnetic 
field. We display clear evidence that this low temperature 
magneto-thermal effect is associated to the existence of low- 
lying fluctuating singlets, but also that the whole picture (T 
behavior of C v and thermally activated spin susceptibility) 
implies contribution of both non magnetic and magnetic ex- 
citations. Comparison with experiments is made. 



Introduction: Spin systems with strong geometri- 
cal frustration exhibit a generic feature: the deviation 
of their spin susceptibility from the Curie law occurs at a 
much lower temperature than for non-frustrated ones and 
is often accompanied by very peculiar low-temperature 
properties [jjj. The two-dimensional spin-i kagome anti- 
ferromagnet (KAFM) is an example of such systems. 

All theoretical and numerical treatments point to a 
disordered Spin-Liquid ground-state with a small 

spin-gap estimated to be of the order of 1 /20 of the cou- 
pling constant. But contrarily to what prevails in 1-D or 
ladders Spin-Liquids, we have shown that the first excita- 
tions of this system have a spin 0. The spectrum of these 
many-particle <Stot = eigen-states forms a continuum 
adjacent to the ground-state, which extends smoothly 
up to energies higher than the spin-gap I n this 

paper we report some thermodynamic consequences of 
this property and focus on the magneto-thermal effects 
which are seemingly a unique signature of this kind of 
Spin-Liquids. 

Numerical method: The i-spins on the Kagome lat- 
tice interact through the Hamiltonian: 



(1) 



where the double sum is limited to first neighbors. The 
temperature T and the magnetic field H are measured 
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in coupling constant units (J = 1,7 = 1). Up to now 
two kinds of methods were available to obtain thermody- 
namic quantities: high temperature series and exact diag- 
onalizations. The first M gives essentially exact results 
down to temperatures of the order of J/2, but unhap- 
pily no direct information on the range of temperature of 
interest in this work. The second allows the exact com- 
putation of the thermodynamics properties at all values 
of T, but for a set of sizes limited to N < 21. For such 
sizes and in spite of the supposed-to-be short range spin- 
spin correlations their might be some uncertainties 
on the low T thermodynamic properties related to the 
discrete nature of the finite size spectra. To deal with 
larger sizes and in order to obtain significant information 
on the effect of an applied field at low T ( T ~ 0.05. .0.2), 
we have devised an hybrid method which takes advantage 
of both previous approaches. An approximate density of 
states in each Irreducible Representation (I. R.) of the 
complete symmetry group of the lattice and Hamilto- 
nian is reconstructed via a maximum entropy procedure. 
It uses exact diagonalization data to fix the edges of the 
density of states and the six first moments of the Hamil- 
tonian in each given subspace. In the N — 36 sample 
there are 264 different I.R. and thus the procedure in- 
volves a very large amount of exact results. Most of the 
low T physics comes from the relative position, weight 
and general shape of the different I.R., which are an es- 
sential input of our method. The knowledge of these 
quantities for each I.R. of SU(2) allows a straightforward 
computation of the effect of any magnetic field. We have 
checked on different frustrated spin problems that this 
approach gives reliable information on the intermediate 
temperature range and on multiple-peak structures j^] . 

Numerical results: High temperature range down 
to T = 0.2: Numerical results are summarized in 
Fig. 1, together with the Pade approximant of the high- 
temperature series of Elstner and Young || . In this range 
of temperature, the spin susceptibility \ shows a shoul- 
der at T ~ 1 and the specific heat C v a high T peak 
around T = 2/3. The Pade approximants cannot be dis- 
tinguished from the N = 18 results down to T = 0.4 
and then merge smoothly with the exact N = 24 low T 
results. We thus infer that above T = 0.4 these different 
data give a very good approximation of the thermody- 
namic limit. It should be noticed that there is only 50% 
of the total entropy in the high T peak above T = 0.2, 
whereas the Neel ordered Heisenberg magnet on the tri- 
angular lattice develops 50% of the total entropy above 
T = 1, and only - 6% below T = 0.2. The low T entropy 
is indeed a distinctive mark of exotic spin liquids J10[ and 
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the KAFM represents an extreme situation. 

From T = 0.2 down to T ~ 0.2A: when cooling down 
the KAFM, one encounters a second energy scale A re- 
lated to the spin-gap pi| . In the range of sizes presently 
studied, the gap is still decreasing fij and A varies from 
Ais = 0.104 to A36 = 0.074 which explains the evolution 
of C v in Fig. 2, 3. From the study of spin-spin correla- 
tions HH and of finite size effects on the spin-gap one 
might expect that the N — 36 sample is not very far from 
the thermodynamic limit. Remembering that 50% of the 
total entropy has to be accounted for below T = 0.2, one 
expects the low T peak in C v to remain for N — > 00. 

The spin-gap clearly explains the thermally activated 
spin susceptibility \- However, C v (in contrast to \) 
does not decrease exponentially in the gap. For the larger 
samples (N = 18 — 36), C v below the peak is well fitted 
by a T 2 law in the range T = 0.3 - 0.6 A (Fig. 2). The 
non exponential behavior of C v is indeed due to the pres- 
ence of singlets in the gap, but, as can be seen in Fig. 3, 
higher spin channels do contribute approximately to one 
half of C v in this range of temperature and the T 2 behav- 
ior appears as a cooperative effect due to different spin 
channels. This fact contradicts an intuitive expectation: 
from our knowledge of fully gapped systems one would 
naively expect an absence of contribution of the triplets 
for temperatures lower than the spin-gap. In fact C v , at 
temperature T, essentially measures the energy fluctua- 
tions around the average energy e(T). In a fully gapped 
system, when T is lower than the spin-gap, the ground- 
state alone is populated and the energy fluctuations are 
thermally activated. In the present case, the thermal 
population at a temperature T is dominated by the sin- 
glet states with an energy ~ e(T): this large number of 
relay levels in the spin-gap and the very large number 
of low lying triplets are at the origin of this surprising 
behavior. 

There has been in the past many speculations about 
the origin of the T 2 behavior in SrCrGaO and in the 
jarosites : simple or nematic spin-waves modes have 
been suggested , . Strictly speaking the presence 
of a spin-gap precludes such an explanation in the spin-4 
case, but the small spin-gap is a proof that the spin- 
i KAFM is relatively near a quantum critical point. 
An alternative explanation can also be featured by fo- 
cussing on the Resonating- Valence-Bound picture of the 
singlets H,||jia): in such a picture the modes are asso- 
ciated to the dimer-dimer long range orientational order 
(resonons of Rokhsar and Kivelson H] ) . Such a picture 
gives a T 2 behavior of C v provided that the dimer-dimer 
structure factor is linear in k at small wave-vectors. This 
would imply an algebraic decrease of the dimer-dimer 
correlations. Regarding the range of sizes available, this 
hypothesis cannot be completely discarded (dimer-dimer 
correlations at distance 6 are of the order of 10~ 2 in the 
spin-i KAFM whereas they are about 10~ 4 in the 
M.S.E. [|l0|). In this last hypothesis, a global picture in- 



volving both singlet and magnetic excitations should be 
worked out. 

Temperatures lower than 0.2A; The C v curves for 
N = 18 and N = 36 are numerically exact in this temper- 
ature range. The entropy of the singlet excitations in the 
spin-gap (~ 15%) is certainly a significant information in 
the thermodynamic limit. But the low-temperature peak 
(dots in Fig. 1, 2 and 3) is unphysical and a consequence 
of the discreteness of the spectra. Large size effects in the 
lowest part of the energy spectrum preclude its extrapo- 
lation for N — ► 00, in particular, they make it impossible 
to decide the exact value of the entropy at T — and the 
very low T behavior of C v in the thermodynamic limit. 

Comparison with experimental results: The 
comparison with experimental results is only indicative, 
insofar as only pseudo-spin- 1 fl5|| , spin-| fl6| |li| and 
spin- 1 J2(],^l| Kagome lattices have been realized and 
studied in the laboratory. 

Pseudo Spin-1: The organic compound m — 
MPYNN.BF4 studied by Wada and coworkers may 
behave at low enough temperatures as a spin-1 KAFM. 
X has been measured in a large range of temperature. 
There are similarities with the present spin-i results: a 
shoulder at T ~ 1, a thermally activated behavior with 
a spin-gap of the order of 7.7 10~ 2 . It might be noticed 
that the experimental gap is large enough compared to 
the spin-i case; half integer and integer spin systems 
might indeed behave differently. 

Spin-Z/ 2 systems: We will concentrate here on 
SrCrGaO which displays the most convincing example 
of a non-conventional low-temperature behavior [fj"7j-[T9[ . 
From the magnetic point of view, it is not clear whether 
SrCrGaO is a good representative of a system of Kagome 
planes (with intervening dilute triangular Heisenberg 
planes) or whether it is better described as a stack of 
2-dimensional slabs of pyrochlore. However, on the qual- 
itative level these two theoretical models display similari- 
ties, absence of long range order, resonating valence-bond 
ground-state, so that a comparison of the experimental 
results for SrCrGaO with our numerical results might 
be justified |||. 

The double peak structure in C v has not be seen in 
these compounds: this is not a surprise at the light of 
present results. In SrCrGaO the Curie- Weiss temper- 
ature is of the order of 50QK, the coupling constant of 
equation (1) is thus of the order of 100K and from results 
shown in Fig. 3 we expect the low-T peak to appear at 
T ~ 0.05 J ~ 5K or below. On the other hand the scale 
for the semi-classical behavior is given by S(S + 1) J: we 
thus predict that the high T peak in SrCrGaO should 
appear around 350-fT, but experimental data only extend 
to AOK. We might notice that in SrCrGaO as in the 
S = i numerical data, the low-T peak in C v accounts for 
half of the total entropy of the spin system ^3|. Endly 
the T 2 behavior of C v is consistent with the experimental 
data. 
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The stronger argument in favor of a parentage between 
the spin-i KAFM and SrCrGaO is the quasi absence of 
sensitivity of C v to applied magnetic fields H as large as 
twice the temperature reported by Ramirez et al. |23| . 

This is a highly surprising result, which cannot be ex- 
plained in any of the following hypothesis: Neel long 
range order, ordinary spin glass, Spin-Liquid with a full 
gap both for magnetic and non magnetic excitations. 

i) Since a Neel-ordered system spontaneously breaks 
the rotational invariance, the ground-state (and the first 
excitations) are superpositions of quasi-degenerate eigen- 
states with different total spin. In such a case, x goes to 
a constant at low T. Though the magnetization per spin 
vanishes, the typical spin is not zero: Sf ot ~ N\{T = 0). 
This gives a H dependent C v . 

ii) On the other hand in a fully gapped Spin-Liquid, 
as soon as the applied field is large enough to match the 
gap, C v is completely washed out, as for example in the 
Spin-Peierls compounds. 

In SrCrGaO a field of 8 to 11 Tesla is large enough to 
excite magnetic excitations (as can be seen in looking at 
the non zero values of \ a t equivalent temperature) and 
nevertheless there is no effect on C v . We show on the 
spin-i example below, that the presence of an impor- 
tant background of non-magnetic excitations is enough 
to drastically decrease the effect of the field on C v . 

The magneto-thermal effect on spin-i KAFM is dis- 
played in Fig. 3 for a field equal to the temperature of 
the peak T p (corresponding approximately to the 6 Tesla 
experiments of Ramirez et al. ) . 

i) The contribution of singlet excitations to C v is dis- 
played on the same graph. In the interesting range of 
temperature, the singlets only account for half of C v . 

ii) The higher spin channels are responsible for the 
missing half of C v . These higher spin channels are in- 
deed fully sensitive to H . If these magnetic levels were 
separated from the singlet ground-sate by a spin-gap void 
of singlet excitations the decrease of the C v peak would 
be 3 times larger than it is in the present situation. Here 
again, the presence of the background of singlets brings 
a great change in the fluctuation spectrum of this system 
and explains its paradoxical low sensitivity to H. 

iii) The effect of a field of the order of T p seems larger in 
the spin-i KAFM than in SrCrGaO: this suggests that 
in this range of T the relative weight of the singlets is 
larger in SrCrGaO than it is for spin-^. We suggest the 
following conjecture: SrCrGaO, as a spin-| system, is 
closer to the quantum critical transition than the spin-^ 
KAFM. Its spin-gap, if not zero is certainly smaller than 
that of the spin-i KAFM, and there may be a collapse 
of the different low-energy scales appearing in the low 
S to t sectors of the spin-i system (see Fig. 3) . In such a 
situation we expect a weaker sensitivity of C v to H. 

An important issue in experimental systems is the 
problem of dilution. We have checked the effect of 
quenched disorder on the spin-^ KAFM by including 1 



or 3 non magnetic sites in a N = 21 sample. The general 
qualitative properties of the KAFM (thermally activated 
X, low-T peak in C v and singlets in the spin-gap) do not 
disappear with such dilutions! This seems an interest- 
ing indication that the thermodynamic features we have 
been describing in this paper are rather robust. 

Conclusion: In this paper we describe the results of a 
numerical investigation of the magneto-thermal proper- 
ties of the spin-i KAFM. We show that the specificities 
of its low energy spectrum induce a very low sensitivity 
of the specific heat to the external magnetic field. Such 
an effect is a very strong signature of a large weight of 
singlet excitations and would not be present neither in a 
Neel ordered system, in a standard Spin-Glass nor in a 
dimerized Spin-Peierls like system. We have verified that 
this effect is qualitatively insensitive to quenched disor- 
der. The qualitative similarities to Ramirez et al results 
on SrCrGaO are probably more than a coincidence but 
a new independent proof advocating for the presence of 
low lying fluctuating singlets in this spin-| compound. 
The precise nature of the lattice (Kagome or pyrochlore 
slab) might not be a central issue. We have preliminary 
results showing that similar Spin-Liquid behavior does 
appear on many other lattices. We now suspect this be- 
havior to be a generic feature of one class of Spin-Liquids. 

Acknowledgments: We thanks A. Ramirez for discus- 
sion and communication of his results a long time before 
publication. We have benefited from very interesting dis- 
cussions with J. P. Boucher, N. Elstner, M. Gingras, S. 
Kivelson and P. Mendels. Computations were performed 
on C98 and T3E at the Institut de Developpement des 
Recherches en Informatique Scientifique of C.N.R.S. un- 
der contracts 984091-980076 and on the T3E at NIC 
Julich. 



[1] A. P. Ramirez, Ann. Rev. Mater. Sci 24, 453 (1994), P. 
Schifler and A. P. Ramirez, Comments Condens. Matter 
Phys. 18, 21 (1996) and references therein. 

[2] V. Elser, Phys. Rev. Lett. 62, 2405 (1989). 

[3] C. Zeng and V. Elser, Phys. Rev. B 51, 8318 (1995). 

[4] J. Chalker and J. Eastmond, Phys. Rev. B 46, 14201 
(1992). 

[5] P. Leung and V. Elser, Phys. Rev. B 47, 5459 (1993). 
[6] P. Lecheminant et al, Phys. Rev. B 56, 2521 (1997). 
[7] C. Waldtmann et al, Eur. Phys. J. B 2, 501 (1998). 
[8] N. Elstner and A. P. Young, Phys. Rev. B 50, 6871 
(1994). 

[9] G. Misguich, Ph.D. thesis, Paris VI, France (1999). 
[10] G. Misguich, B. Bernu, C. Lhuillier, and C. Waldtmann, 
Phys. Rev. Lett. 81, 1098 (1998) and Phys. Rev. B 60, 
1064 (1999). 

[11] P. Chandra and P. Coleman, Phys. Rev. Lett. 66, 100 
(1991). 



3 



[12] I. Ritchey, P. Chandra, and P. Coleman, Phys. Rev. B 

47, 15342 (1993). 
[13] F. Mila, Phys. Rev. Lett. 81, 2356 (1998). 
[14] D. Rokhsar and S. Kivelson, Phys. Rev. Lett. 61, 2376 

(1988). 

[15] N. Wada et al, J. Phys. Soc. Jpn. 66, 961 (1997). 
[16] X. Obradors et al, Solid State Commun. 65, 189 (1988). 
[17] A. Ramirez et al, Phys. Rev. Lett. 64, 2070 (1992). 
[18] Y. Uemura et al, Phys. Rev. Lett. 73, 3306 (1994). 
[19] S.-H. Lee et al, Europhys. Lett 35, 127 (1996). 
[20] A. Keren et al, Phys. Rev. B 53, 6451 (1996). 
[21] A. S. Wills et al, Europhys. Lett 42, 325 (1998). 
[22] B. Canals and C. Lacroix, Phys. Rev. Lett. 80, 2933 
(1998). 

[23] A. P. Ramirez, B. Hessen and M. Winklemann (to appear 
in Phys. Rev. Lett.) 

[24] For each size, this energy scale A is obtained from a ther- 
modynamic approach of the energy per spin e vs magne- 
tization m — Stat/N: e(m) ~ e(0) + Am + am 2 at small 
m. 



0.15 



X 



Z 
> 

o 




0.05 - 



FIG. 1. Specific heat C v , spin susceptibility \ and entropy 
S of the N — 18 Kagome antiferromagnet versus temperature 
T. Dashed line: C v /Nks ; continuous line: x x 2; short-dashed 
line: 5 , /fcsln2. The light dashes represent the Pade approxi- 
mants to C v and \ °f Ref- [§ : they cannot be distinguished 
from the N = 18 and N = 24 results down to T = 0.4. For 
T ^ 0.03 the numerical results are exact but cannot be infor- 
mative on the thermodynamic limit, because of the finite size 
effect. 
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FIG. 2. Low-temperature specific heat C v (full lines for 
T > 0.3A, dotted lines below ) versus the square of the tem- 
perature (JV = 18 and 36 samples). C v is well fitted by a T 2 
law (light line) in the range T = 0.3 — 0.6A and deviates from 
higher or lower power laws. 
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FIG. 3. Contribution of the singlet (Stot = 0) states (dot- 
ted lines) to the total specific heat (full lines) and effect of 
a magnetic field H equivalent to the 6 Tesla experiment of 
Ramirez et al (dashed lines). The arrows indicate the 
value of A for each sample. 
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